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We analyze the entanglement of SU(2)-invariant density matrices of two spins Si, S2 using the Peres- 
Horodecki criterion. Such density matrices arise from thermal equilibrium states of isotropic spin 
systems. The partial transpose of such a state has the same multiplet structure and degeneracies as 
the original matrix with the eigenvalue of largest multiplicity being non-negative. The case Si = S, 
S2 = 1/2 can be solved completely and is discussed in detail with respect to isotropic Heisenberg spin 
models. Moreover, in this case the Peres-Horodecki ciriterion turns out to be a sufficient condition 
for non-separability. We also characterize SU(2)-invariant states of two spins of length 1. 



I. INTRODUCTION 

Entanglement is one of the most intriguing properties 
of quantum physics [1,2] and the key ingredient of the 
emerging field of quantum information theory and pro- 
cessing [3]. Very recently substantial interest has devel- 
oped concerning the question of quantum entanglement 
in equilibrium states of quantum spin systems [7-22] as 
often studied in condensed matter physics and statistical 
mechanics. More specifically, a typical question arises: 
Tracing out from a many-body system all degrees of free- 
dom except for, say, two spins, is this reduced density 
matrix separable or not? 

In the present work we study SU(2)-invariant density 
matrices of two spins. These states are are defined to 
be invariant under all uniform rotations U\®U2 of both 
spins S\ and S2, where U a = exjp(irfS a ), a G {1,2}, are 
transformations corresponding to the same set of real pa- 
rameters fj in the representation of SU(2) appropriate for 
the spin lengths S\ and S% {h= 1). In other words, such 
states p commute with all components of the total spin 
J = S%+ S%- SU(2)-invariant density matrices arise from 
thermal equilibrium states of spin systems with an rota- 
tionally invariant Hamiltonian by tracing out all degrees 
of freedom but those two spins [23] . Our results general- 
ize previous work on such systems [7,8,13,14] to the case 
of higher spins, and we discuss our findings with respect 
to generic antiferromagnetic or ferrimagnetic Heisenberg 
spin lattice models. These considerations lead to the nat- 
ural conclusion that pairwise quantum entanglement in 
equilibrium states is strongest in systems with small spin 
length and low spatial dimension. Moreover, an insepa- 
rable equilibrium reduced two-spin state can usually only 
be achieved for neighboring spins, but not for more dis- 
tant lattice sites. Viewed like this, equilibrium states of 
such systems do not appear to be a particularly strong 
source of pairwise quantum entanglement. 

To investigate the separability of this type of density 
matrices in the case of higher spins we shall make use 
of the Peres-Horodecki criterion [24,25]. This criterion 
states that a separable density matrix has necessarily a 
positive partial transpose [24]. Moreover a positive par- 



tial transpose is also sufficient [25] for the separability 
of a given density matrix in the case of two qubits, and 
in the case of a qubit and a qutrit (i.e. the case of two 
spins of length 1/2, and of a spin 1/2 and a spin 1, re- 
spectively). For larger dimensions of the parties (spins) 
involved non-separable states with positive partial trans- 
pose exist [26]. To apply the Peres-Horodecki criterion 
one can perform the partial transposition with respect 
to either subsystem (spin), since both resulting matrices 
have the same spectrum. For definiteness we will consider 
in the following the partial transpose p T2 with respect to 
S 2 . 

SU(2)-invariant two-spin states where mentioned 
briefly already by Vollbrccht and Werner [27], where it 
was pointed out that the case S± = S2 = 1 corresponds 
to states invariant under O ®0 with O being a 0(3) ro- 
tation. In this work we will also give explicit criteria in 
terms of spin correlators for such states to have a positive 
partial transpose [28] . 

II. SU(2)-INVARIANT STATES AND THEIR 
PARTIAL TRANSPOSE 

Let us start with some general considerations. Since 
an SU(2)-invariant state commutes with all components 
of J, it acts, according to Schur's Lemma, as a scalar on 
each irreducible representation (multiplet) of J. There- 
fore p has the general from 

p= J2 WTi E \J,J z )oo(J,J z \, (i) 

J=\S x -S a \ J*=-J 

where the constants A(J) fulfill A(J) > 0, J2j A ( J ) = 
1. Here \ J,J z )o denotes a state of total spin J and z- 
component J z . 

Now let O be a general operator acting on a bipartite 
system. If O is transformed by U\ <g> U2 , T2 transforms 
covariantly under U\ <E> U% ■, 

({u 1 ®u 2 )o(u+ ®u+)) T2 

= {u l ®u* 2 )o T2 (u+®{u;) + ) . (2) 
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Here U\, U2 are general unitary transformations acting 
on the subsystems and do not necessarily represent SU(2) 
transformations. The relation (2) can be derived readily 
by writing O in the form Oi 1 i 2 j 1 j 2 where the subscripts 
at the indices refer to the different subsystems [27,29]. In 
particular, it follows that if O is invariant under U\ ® U2, 
T2 is invariant under U\ ® U%- 

Now let Ui, U2 represent again SU(2) transformations. 
In the standard representations of their generators the x- 
and z-components are given by real matrices while the 
matrices for the y-components are imaginary. Thus, a 
complex conjugation of C/ 2 = exp(irfS 2 ), is equivalent 
to changing the sign of Sf, Sf . Therefore, p T2 com- 
mutes with the operators K defined by K x = Sf — Sf , 
K v = Sf + S 2 , K z = Sf - Sf , and these operators also 
furnish a representation of su(2), [if", if' 3 ] = ie a ^ 1 K 1 
(using standard notation). In the basis of tensor prod- 
uct states |Sf,Sf) of Sf, Sf, p T2 is block-diagonal with 
respect given values of K z . In particular, | ± Si,+S 2 ) 
are eigenstates of p T2 with the degenerate eigenvalue 
(±Si, =fS 2 |/0|±Si, =fS 2 ) > 0. Now it follows from elemen- 
tary representation theory that p 7 " 2 has actually an SU(2) 
multiplet structure with respect to the operators K. The 
multiplets are labeled by the value of K 2 = K(K + 1) 
with I Si - S 2 | < K < (Si + S 2 ). On these multiplets 
p T2 acts as a constant. As seen above, the eigenvalue 
corresponding to the largest if- multiplet is always non- 
negative. 



III. THE CASE £2 = 1/2 

Let us now consider system consisting of a spin Si of 
arbitrary length S and a spin S 2 of length 1/2. Here a 
general SU(2)-invariant density matrix has the form 



S-1/2 



E l^-i J z )oo(S-i,J 

J* = -S+l/2 
S+l/2 



+ 



2S 

1 - F 
2S + 2 



E \S+\,J z )m{S + \,J z \. 
-s-1/2 



(3) 



The quantity F £ [0, 1] is in thermal equilibrium a func- 
tion of temperature and, in the case of p being a reduced 
density matrix of a larger system, it contains information 
about the entire system which has been traced out ex- 
cept for the spins Si and S 2 . By expressing F in terms 
of the projector onto the J = S — 1/2 multiplet one finds 
F = (S - 2(SiS 2 ))/(2S + 1), where (•} denotes the ex- 
pectation value with respect to p. Thus, p is completely 
determined by the correlator (SiS 2 ). 

In order to perform a partial transposition on p it is 
convenient to express it in a basis of tensor product eigen- 
states |S z ,±l/2) of Sf and Sf. Using the well-known 



Clcbsch-Gordan coefficients for coupling a spin S to a 
spin 1/2, the nonvanishing matrix elements are given by 



<S z ,±l/2|HS z ,±l/2> 

1 / (StS z )F (S ± S z + 1)(1 - F) 
~ 2S+ 1 ^ 2S + 2S + 2 

(S z , l/2\p\S z + 1,-1/2) 

= y/(S-S z )(S + S^ + l) 

2S + 1 



F 1 — F 

"2S + 2S + 2 



(4) 



(5) 



The partial transpose p T2 is diagonal on the sub- 
space spanned by S, — 1/2), | — S, +1/2) lying in 
the largest K multiplet with the eigenvalue Ai := 
(S, -1/2HS, -1/2) = (-S, l/2|p| - S, 1/2) where 



Ai = 



1 



2S+ 1 



F 



1 - F 
2S + 2 



(6) 



On the remaining Hilbcrt space the partial transpose 
is block-diagonal, where the blocks act on subspaces 
spanned by the basis vectors |S z ,-l/2), |S Z + 1,1/2) 
and have the from 



<S z ,-l/2|p|SV 
(S*,l/2|p|S z + l, 



1/2) <S*,1/2|HS* + 1,-1/2) 
-1/2) (S* + l,l/2|p|S* + 1,1/2) 



The eigenvalues of these submatrices are given by Ai and 
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2S + 1 2S 



(7) 



These eigenvalues do not depend on S z . Therefore Ai 
and A 2 occur with the multiplicities 2S + 2 and 2S, re- 
spectively, in accordance with the above general results. 
Moreover, Ai is always positive, while A 2 becomes nega- 
tive for F > 2S/(2S + 1), or, equivalently, 



(S!S 2 )<--. 



(8) 



Thus, our state has non-positive partial transpose if and 
only if the correlator (SiS 2 ) is negative and larger in 
modulus than S/2. This is the maximum value |(SiS 2 )| 
achieve in a separable state. This intuitively very rea- 
sonable criterion includes earlier results by Wang and 
Zanardi [14] who investigated the case Si = S 2 = 1/2 
by evaluating the entanglement of formation [30] using 
Wootters' concurrence [31]. Unfortunately this is not 
a viable route for Si > 1/2 since Wootters' construc- 
tion appears to be restricted to the case of two qubits. 
Moreover, with increasing S, the states with non-positive 
partial transpose have increasing weight in the smaller 
multiplet J = S — 1/2, approaching unity for S — > 00. 

If a given state has a negative partial transpose it is 
necessarily entangled. Moreover, in the case of SU(2)- 
invariant states with S 2 = 1/2 studied in this section, a 
positive partial transpose, i.e. (SiS 2 ) > —S/2, is also a 
sufficient criterion for separable states. We prove this fact 
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by explicitly constructing a decomposition consisting of 
projectors on porduct states. If a given state has a posi- 
tive partial transpose we can write (/Si 6*2) = (S/2) cos(7) 
with some real angle 7. Now let |0) denote a spin- 
coherent state [32,33] of Si pointing in some arbitrary 
direction, and I7) denote a spin-coherent state of S2 with 
its polarization direction forming the angle 7 with the 
polarization direction of Si. In the pure product state 
|0) <8> I7) we have by construction (S1S2) = (S/2) cos(7), 
and the value of this correlator is invariant under all uni- 
form rotations of both spins. Now consider 



p = N(S) J d 3 V [(U 1 (f 1 )<E>U 2 (v)) 

(|o>® |7><o|® {i\){Ui(ii)®umy 



(9) 



where the integration goes over all simultaneous rotations 
parametrized as U a — exp(iffS a ), a <E {1, 2}, and N(S) is 
a normalization constant. Thus the state (9) is a separa- 
ble state which is obviously invariant under simultaneous 
rotations of both spins and fulfills (S1S2) = (S/2) cos(7). 
Since such an SU(2)-invariant state is uniquely deter- 
mined by this correlator, we have constructed a decom- 
position of the original state in terms of projectors on 
product states, which completes the proof. 

Let us now discuss the above result with respect to 
isotropic Heisenberg lattice spin models as studied inten- 
sively in condensed matter physics and statistical me- 
chanics. Previous studies have concentrated on one- 
dimensional systems. This has on the one hand the 
practical reason that for such systems the body of exact 
results concerning correlations is largest. On the other 
hand this is due to the fact that quantum correlations can 
generically be expected to become weaker with increasing 
spatial dimension, i.e. with increasing number of neigh- 
bors to each spin. Therefore one-dimensional systems 
are the most attractive to look for equilibrium quantum 
entanglement. 

Since quantum correlations like (S1S2) can generally 
be expected to decay with increasing temperature, the 
criterion (8) defines implicitly a threshold temperature 
for the occurrence of a non-positive partial transpose, 
provided the inequality (8) is fulfilled in the ground state 
at T = [14]. This can only be the case in antiferromag- 
nctic or, for Si > S2 = 1/2, fcrrimagnetic systems. In 
particular, in an antiferromagnetic spin 1/2 chain corre- 
lations arc generically of the from 



{S m S m+n ) = (-iy i x(n) , 



(10) 



where m denotes some lattice site in the translation- 
ally invariant chain, n is the number of lattice sites 
between the spins considered, and x is a positive and 
monotonously decaying function. The alternating sign 
resembles Neel ordering as it is present in the ground 
state of a classical antiferromagnet. It follows that the 



reduced equilibrium two-spin density matrix can only be 
entangled if the spins involved reside on different sub- 
lattices (corresponding to odd n). This intuitively clear 
finding also holds for generic antiferromagnets or ferri- 
magnets on bipartite lattices in higher spatial dimension. 

The spin 1/2 Heisenberg chain with antiferromag- 
netic exchange between nearest neighbors is described by 
the Hamiltonian Ti — J2 m SmS m +i Here the correlator 
(S m S m+ i). is equal to the ground state energy per spin 
and given by In 2 - 1/4 « -0.443 < -0.25 [34]. Thus, 
the criterion (8) is fulfilled [14]. However, for larger dis- 
tances between the spins, n 6 {3, 5, . . .}, numerical data 
[35] shows that the inequality (8) is violated and, accord- 
ing to the Peres-Horodecki criterion, the corresponding 
reduced density matrix is separable. The same state- 
ments apply to another typical antiferromagnetic spin 
1/2 wave function, the so-called Gutzwiller wave func- 
tion [36], where {S m S m+n ) can be evaluated analytically 
for all n [37]. 

In summary, one-dimensional antiferromagnetic iso- 
tropic Heisenberg models of spins 1/2 do generically 
not appear to be in thermal equilibrium a particularly 
strong source of entanglement, since usually only the re- 
duced density matrices of neigboring spins are insepara- 
ble, while all others are non-entangled. This result might 
appear somewhat surprising since such systems are usu- 
ally considered to have particularly strong spin corre- 
lations because of the small spin length as well as the 
low spatial dimension [33]. However, as seen here, these 
strong quantum correlations do in general not translate 
to long-ranged entanglement in equilibrium reduced den- 
sity matrices. 

As pointed out already, this finding cannot be expected 
to change in the case of higher spatial dimension. More- 
over, ferrimagnetic systems involving spins Si > 1/2 will 
generically have the same properties, since we have the 
same type of criterion (8) for the partial transpose being 
non-positive. This criterion requires sufficiently strong 
quantum fluctuations which are generically reduced with 
increasing spin length. To illustrate these trends let us 
rewrite the criterion (8) in the form (S1S2} /(S1S2) < — 1- 
For a ferrimagnetic chain consisting of alternating spins 
Si = 1, S2 = 1/2 a numerical estimate for the l.h.s. of 
this inequality in the case of neighboring spins is given 
by [38] (SiS 2 )/(SiS 2 ) = -1.455. For a two-dimensional 
spin- 1/2 antiferromagnet on the square lattice one finds 
[39] (SiS 2 )/(SiS 2 ) = -1.344. Both values are larger 
than the result for the spin-1/2 chain as discussed above, 
(SiS 2 )/(SiS2) = —1.773, indicating the suppression of 
pairwise entanglement with increasing spatial dimension 
and lengths of spins involved. 
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IV. THE CASE S 2 > 1 

We now turn to the case 5i =: S > 1, 5 2 = 1. Here 
the general SU(2)-invariant density matrix reads 

s-i 

]T |S-l,J*>oo<S-l,J*| 

J*=-S+l 

s 



p = 



G 
25-1 



if 



+ ^5+T S |5,J z >oo<S,J*| 



+ 



j*=-s 
1-G-H 



s+i 



5 + 1, J*) 00 <S + 1, J*|. (11) 



Expressing the quantities G and if in terms of projectors 
onto the multiplets on the total spin Je {5— 1,5, 5+1} 
one finds 

° = S(2S + 1) t S ~ {S ~ 1)(S?1 ^ 2> + (12) 
H = l- j^—^ (0i&) + ((SiS 2 f)) . (13) 

Therefore, p is completely determined by the correlators 
(S1S2) and ((5i5 2 ) 2 ) and has three different eigenvalues 
with degeneracies 25— 1, 25+1, and 25 + 3, correspond- 
ing to multiplets of the total spin J. The three different 
eigenvalues of p T2 can be found in the subspaces spanned 
by (\S* + 1,1>,|S*,0),|S* - 1 - 1}), \S Z \ ± 5 having 
K z = S z . With respect to this basis p T2 reads 

a(S z ) 5(S Z ) e(S z ) 

S(S Z ) (3(S Z ) r,(S z ) J . (14) 

s(S z ) V (S Z ) j(S z ) 

with 



a(S z ) = (S z + l,l|p|5 2 + 1,1) 
0(S z ) = (S z ,O\p\S z ,O) 
7 (5^) = (5 z -l-l|p|5 2 -l,-l) 
S(S') = {S',l\p\S' + l,0) 
e(S') = (S' + l,-l\p\S'-l,l) 
V (S z ) = (S z ,-l\p\S z -l,0) 



(15) 
(16) 
(17) 
(18) 
(19) 
(20) 



Unfortunately, the evaluation of the matrix elements in 
(14) for general values of 5 and S z turns out to be ex- 
tremely tedious because the form of the Clebsch-Gordan 
coefficients is drastically more complicated than in the 
previous case 5 2 = 1/2. For simplicity we therefore con- 
centrate on the case 5=1 where the only block of the 
form (14) corresponds to S z = with a — 7 and S = r/. 
In this case we have the eigenvalues 





1 


3 


Ml = 


30 " 


"To 




1 




^2 = 


6 ~ 






1 




fJ-3 = 


3 ~ 


3 



(21) 
(22) 
(23) 



The eigenvalue p\ is always positive and corresponds to 
the largest if-multiplet. Thus, p T2 has negative eigen- 
values if and only ifG>l/3or_ff>l/2. Expressed in 
terms of correlators these conditions read 



2 < {{S1S2?) 

1 > (5!5 2 ) + ((5!5 2 ) 2 ) 



(24) 
(25) 



Similarly as in the previous case S\ = 5, 5 2 = 1/2, 
the two-spin density matrix has a non-positive partial 
transpose only if the weight of the smaller multiplets is 
sufficiently large. We note that, differently from the case 
Si = 5, 5 2 = 1/2, an analogous proof for the sufficiency 
of the Peres- Horodecki criterion cannot be given for 5 2 > 
1 since such states arc determined by more than just a 
single correlator. 



V. CONCLUSIONS 

We have analyzed the pairwise quantum entanglement 
in SU(2)-invariant quantum spin systems by studying 
their behavior under partial transposition. As a general 
result, the partial transpose of an SU(2)-invariant two- 
spin state has the same multiplet structure and degenera- 
cies as the original matrix with eigenvalue of largest mul- 
tiplicity being non-negative. SU(2)-invariant density ma- 
trices arise from thermal equilibrium states of isotropic 
spin systems in sufficiently low spatial dimension. The 
case Si = 5, 5 2 = 1/2 can be solved completely and is 
discussed in detail with respect to well-known Heisenberg 
spin models. Moreover, in this case the Peres-Horodecki 
ciriterion turns out to be a sufficient condition for non- 
separability. As a general trend, low spatial dimension 
and small lengths of the spins involved tend to facilitate 
the occurrence of inseparable equilibrium states. How- 
ever, inseparability occurs typically only between neigh- 
boring spins in such spin lattice systems. In this sense, 
isotropic Heisenberg spin models do not appear to be 
a particularly strong source of quantum entanglement, 
at least as far as their equilibrium properties are con- 
cerned. Finally we have also characterized the properties 
of SU(2)-invariant states of two spins of length 1 under 
partial transposition. 
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